Possible Emission of Cosmic X— and 7— rays by Unstable Particles at Late Times 
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Not all astrophysical mechanisms of the emission of electromagnetic radiation including X- and 
7- rays coming from the space are clear. We find that charged unstable particles as well as neutral 
unstable particles with non-zero magnetic moment which live sufficiently long may emit electro- 
magnetic radiation. This new mechanism is connected with the properties of unstable particles at 
the post exponential time region. Analyzing the transition time region between exponential and 
non-exponential form of the survival amplitude it is found that the instantaneous energy of the 
unstable particle can take very large values, much larger than the energy of this state for times from 
the exponential time region. Basing on the results obtained for the model considered, it is shown 
that this purely quantum mechanical effect may be responsible for causing unstable particles to emit 
electromagnetic-, X— or 7-rays at some time intervals from the transition time regions. 
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Typical physical processes in which cosmic microwave 
and other electromagnetic radiation, X—, or 7-rays are 
generated have purely electromagnetic nature (an accel- 
eration of charged particles, inverse Compoton scatter- 
ing, etc.), or have the nature of nuclear and particle 
physics reactions (particle— antiparticle annihilation, nu- 
clear fusion and fission, nuclear or particle decay). The 
knowledge of these processes is not sufficient for explain- 
ing all mechanisms driving the emission from some galac- 
tic and extragalactic X- or 7-rays sources, e.g. the 
mechanism that generates 7-ray emission of the so-called 
"Fermi bubbles" remains controversial, the mechanisms 
which drive the high energy emission from blazars is still 
poorly understood, etc. (see eg. [H-Q). Astrophysical 
processes are the source of not only electromagnetic, X- 
or 7-rays but also a huge number of elementary particles 
including unstable particles of very high energies (see eg. 
Q). The numbers of created unstable particles during 
these processes are so large that some of them can sur- 
vive up to times t at which the survival probability de- 
pending on t transforms from the exponential form into 
the inverse power-like form. It appears that at this time 
region a new quantum effect is observed: A very rapid 
fluctuations of the instantaneous energy of unstable par- 
ticles take place. These fluctuations of the instantaneous 
energy should manifest itself as fluctuations of the veloc- 
ity of the particle. We show that this effect may cause 
unstable particles to emit electromagnetic radiation of a 
very wide spectrum: from radio- up to ultra-high fre- 
quencies v including X-rays and 7-rays. 

Let us start from a brief introduction into the problem 
of the late time behavior of unstable states. If \<f>) is an 
initial unstable state then the survival probability, V(t), 
equals V(t) = \a(t)\ 2 . where a(t) is the survival ampli- 
tude, a(t) = t), and |0;t) = e~ ltH |</>), H is the to- 
tal Hamiltonian of the system under considerations. The 
spectrum, a(H), of H is assumed to be bounded from be- 
low, a{H) = [E m i n , 00) and E m i n > -co. Searching fol- 
iate time properties of unstable states one usually uses 
the integral representation of a(t) as the Fourier trans- 



form of the energy distribution function, ui(E), 



a(t) 



uj{E)e- itE dE, 



(1) 

with w(E) > and oj(E) = for E < E min 
In the case of quasi-stationary (metastable) states it is 
convenient to express a(t) in the following form 111 Il2j. 
a(t) = a exp (t) + ai t (t), where a exp (t) is the exponential 
part of a(i), that is a exp (t) = N exp [-it{Efy - f-T^)], 
(E^ is the energy of the system in the state \4>) mea- 
sured at the canonical decay times, i.e. when V^,{t) has 
the exponential form, E^ is the decay width, N is the 
normalization constant), and au(t) is the late time non- 
exponential part of a(t). 

From the literature it is known that the characteristic 
feature of survival probabilities V{t) is the presence of 
sharp and frequent fluctuations at the transition times 
region, when contributions from |a e ™(t)| 2 and \ait{t)\ 2 
into V{t) are comparable (see, eg. [5|, LZHlfJ), and that 
the amplitude au(t) and thus the probability V{i) ex- 
hibits inverse power-law behavior at the late time region 
for times t much later than the crossover time T. (This 
effect was confirmed experimentally not long ago |14j). 
The crossover time T can be found by solving the follow- 
ing equation, \a exp (t)\ 2 = |a ;t (i)| 2 . In general T > r , 
where = 1/7^, is the live-time of 4>. Formulae for T 
depend on the model considered (i.e. on u>(E)) in general 
(see, eg. EMI ) ■ 

It is commonly known that the information about the 
decay rate, J 1 ^, of the unstable state \<j>) under considera- 
tions can be extracted from the survival amplitude a(t). 
In general not only but also the instantaneous energy 
£d>(i) of an unstable state \<j>) can be calculated using a(t) 



lAuM- I n the considered case, £$(t) can be found using 



the effective Hamiltonian, h^t), governing the time evo- 
lution in an one-dimensional subspace of states spanned 
by vector 
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i da(t) _ (<p\H\(f>;t) 



a{t) dt 



{4>W) 



(2) 
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The instantaneous energy £,f,(t) of the system in the state 
\4>) is the real part of h<p(t), £^,(t) = K (h^(t)). The imag- 
inary part of h$(t) defines the instantaneous decay rate 

r^t), r^ = rj(t) = -2^(h <l> (t)). 

There is £$(t) = E% and r$(t) = E^ at the canonical 
decay times and at asymptotically late times, 



£(f>(t) — E rn i n "h ,2 ^4 



£1 
t 



£3 



(for i > T), (3) 
(for OT), (4) 



where c, = c*, i = 1, 2 , . . ., so lim^oo £$(t) 
limt^oo r (t) =0 [Hill El- Results © and Q arc 
rigorous. The basic physical factor forcing the amplitude 
a(t) to exhibit inverse power law behavior at t 3> T is 
the boundedness from below of cr(H). This means that if 
this condition is satisfied and /_ °° uj(E) dE < oo, then 
all properties of a(t), including the form of the time- 
dependence at t ^> T, are the mathematical consequence 
of them both. The same applies by © to the properties 
of h^(t) and concerns the asymptotic form of h^t) and 
thus of £<f,(t) and E^t) at t > T. 

The sharp and frequent of fluctuations of V$(t) at the 
transition times region are a consequence of a similar be- 
havior of real and imaginary parts of the amplitude a(t) 
at this time region. Therefore the derivatives of a(t) may 
reach extremely large negative and positive values for 
some times from the transition time region and the mod- 
ulus of these derivatives is much larger than the modulus 
of a(t), which is very small for these times. This means 
that at this time region the real part of h^(t) which is 
expressed by the relation ([2]), ic. by a large derivative of 
a(t) divided by a very small a(t), can reach values much 
larger than the energy E® of the unstable state measured 
at the canonical decay times. Using relations ([1]), (TJ|) and 
assuming the form of uj(E) and performing all necessary 
calculations numerically one can see how this mechanism 
work. A typical behavior of the instantaneous energy 
£ <p{t) at the transition time region is presented in Figs (p} 
and (J2J) . In these figures the calculations were performed 
for E min = and the Breit-Wigner energy distribution 

function, u(E) = u BW (E) d = ^ 6(E) {E _ E oyl (r o /2 y , 
where 0(E) is the unit step function. 

Note that from the point of view of a frame of refer- 
ence in which the time evolution of the unstable system 
was calculated the Rothe experiment as well as the pic- 
ture presented in Figs ([I]) , ([2]) refer to the rest coordinate 
system of the unstable system considered. Astrophysical 
sources of unstable particles emit them with relativistic 
or ultra-relativistic velocities in relation to an external 
observer so many of these particles move in space with 
ultra high energies. The question is what effects can be 
observed by an external observer when the unstable par- 
ticle, say </>, which survived up to the transition times re- 
gion, t ~ T, or longer is moving with a relativistic veloc- 
ity in relation to this observer. The distance d from the 
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Figure 1. (a) The instantaneous energy £$(t) in the tran- 
sitions time region: The case E^/T^ = 20. Axes: y = 
£^(t)/E^, x = t/T<i,. The dashed line denotes the straight 
line y — 1. (b) Enlarged part of (a): The highest maximum 
of £ ( / > (t)/E3, in the transition times region. 




Figure 2. The same as in Fig (1) for FA/r? = 10 



source reached by this particle is of order d <~ dy, where 
d T = v*-r, r = 7i Tand 7L = >y L (v+) = (^T^)-\ 
j3 = vv/c, is the velocity of the particle <j>. (For sim- 
plicity we assume that there is a frame of reference com- 
mon for the source and observer both and that they do 
not move with respect to this frame of reference). The 
last expression in relation (J^J) explains why effects of type 
©, (@| and those one can see in Figs (JXJ) , (|2J) arc possible. 
In the case of moving particles created in astrophysical 
processes one should consider the effect shown in Figs 
(Q]), pi together with the fact that the particle gains ex- 
tremely huge kinetic energy, W®, which have to be con- 
served. There is = mj c 2 jl, where to° is the rest 
mass of the particle <fi. We have m" c 2 = E® at canonical 
decay times and thus = E 1 } at these times. At 
this time region i?° = £$(£) but at times t ^> t^, t ~ T 
we have £$(t) ^ E^. So at transition times region the 
relation £<f,(t) = ?7i°(i)c 2 defines the instantaneous rest 
mass 7n°(t) and the kinetic energy of the particle <f> at 
t <~ T equals 

W*{t) = ml(t)c 2 7£ (t) = £^(t) lL (t), (fart~T). (5) 

Here jL(t) = jL(v^(t)), v^(t) is the velocity of the par- 
ticle 4> a t the instant t. Of course the kinetic energies 
W^(t) of <j> have to be the same at the canonical de- 
cay times region and at times t ~ T: = W^(t), that 
is there should be 



= £<p(t)j L (t) = const. 



(6) 
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From relation ([6]) one can infer that this is possible only 
when the changes of £s(t) at times t ~ T are balanced 
with suitable changes of 7i(t) (i.e. of the velocity v^(t) 
of the considered particle). So, in the case of moving 
unstable particles, an external observer should detect 
rapid fluctuations (changes) of their velocities at dis- 
tances d ~ g?t from their source. These fluctuations of 
the velocities mean for the observer that the particles are 
moving with a nonzero acceleration in this space region, 
0. So we can expect that this observer will register 
electromagnetic radiation emitted by charged unstable 
particles, which survived up to times t ~ T, i.e. which 
reached distances d ~ dr from the source (see Fig |3])). 
This follows from the Larmor formula and its rclativis- 
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Figure 3. Time regions: (i) Canonical decay, (ii) Transition, 
(iii) Asymptotically late. Wf = W*(t»), (i = 1,2,3) and 
Wf{ti) is the energy of moving relativistic particle <j>, t\ <C 
t 2 < t 3 , ti ~ T4,, ti ~ T, t 3 > T and W*(ti) = m° c 2 7l = 
E&Il, v is the frequency of the emitted electromagnetic rays. 

tic generalization, which state that the total radiation 
power, P, from the considered charged particle is pro- 
portional to (u^) 2 , P oc (u*) 2 (see e.g. [lj§) and 7^ 
implies that there must be P ^ 0. The same conclusion 
also concerns neutral unstable particles with non-zero 
magnetic moment [l6, 17 1. One should expect that the 
spectrum of this radiation will be very wide: From high 
radio frequencies, through X-rays up to high energy 7- 
rays depending on the scale of the fluctuations of the 
instantaneous energy £s(t) in this space region. 

Within the model defined by lubw(E) the cross-over 
time T can be found using the following approximate 
relation valid for E°J T9 » 1, [ll|: 



T 



2 In 



Mrf) 



(T) 

131 ] one has 
). Consider- 
ing a meson //-"- as an example and using ([7]) one finds 
T = ~ 170 Tjj. (The formula (11) from []J] gives 
T M ~ 202 r M ). The distance dr from the source reached 



whereas for the model considered in 
T/ts ~ 5 In (E°/r°) (see (11) in 




170 t m 
from 



by muon, which survived up to the time T M ~ 
depends on its kinetic energy and equals 
dr„ ^ 10 6 [m] if W = 10 9 [eV], up to d Th ^ 0,027 
[pc] if IF' 1 = 10 18 [cV]. Similarly, for 7r-mesons one 
obtains: dr, 8.5 x 10 3 [m] if W 7 ' = 10 9 [cV] and 
d T ~ 8.6 x'lO 12 [ml if W" = 10 18 [eVl. For the neu- 



tron one finds: d Tn ^ 16, 23 [au] if W n = 10 9 [eV] and 
d Tn ~ 2228 [kpc] if W n = 10 18 [cV]. 

Let us now analyze Fig (1) in more details. Coor- 
dinates of the highest maximum in Fig ([T]) are equal: 
(xmx,Dmx) = (21.60,10.27). Coordinates of points of 
the intersection of this maksimum with the straight line 
y = 1 are equal: (#1,2/1) = (21.58,1.0) and (#2,2/2) = 
(21.62, 1.0). From these coordinates one can extract the 
change Av* = v^(t mx ) - u^(ti) of the velocity of the 
considered particle and the time interval At = t mx — 1\ at 
which this change occurred. Indeed, using ^ one finds 

£(p (j^mx ) 



7i(*i) 



£<j>\k\ 



Ihitrax)- 



(8) 



There is £s(ti) = £7? in the considered case. This means 
that we can replace 7 l(<i) by 7 l measured at the canon- 
ical decay times and then taking the value of the ra- 
tio £<f,(t mx )/E® from Fig ([T]) we can use ([5]) to calculate 
7i(*mx)- Next having 7 z,(ti) = 7l and 7z(i mx ) it is 
easy to find Av^ = v^(t mx ) — v^(t\). Now using rel- 
ativistic generalization of the Larmor formula (see eg. 
formula (14.43) in [l6j]) one can estimate the energy P 
of the electromagnetic radiation emitted in unit of time 
by an unstable charged relativistic particle <fi during the 
time interval At: That is one can find Av^ /At and thus 
P oc (At// At) 2 . This procedure, formulae ©, © and 
parameters describing the highest maximum in Fig. ([1]) 
lead to the following (simplified, very conservative) esti- 
mations of the energies of the electromagnetic radiation 
emitted by ultra relativistic muon at the transition times 
region (in a distance d ~ d? from the source): P ~ 4.6 
[eV/s]. Analogously coordinates of the highest maximum 
in Fig ([2]) are equal: (x mx , y mx ) = (18.69,37.68) and 
coordinates of points of the intersection of this maksi- 
mum with the line y = 1 are: (xi, yi) = (18.67, 1.0) and 
(#2,2/2) = (18.72,1.0). This leads to the following esti- 
mation: P ~ 0.84 [keV/s]. Similar estimations of P can 
be found for neutral ultra-rclativistic unstable particles 
with non-zero magnetic moment. 

The question is where the above described effect may 
be observed. Astrophysical and cosmological processes 
in which extremely huge numbers of unstable particles 
are created seem to be a possibility for the above dis- 
cussed effect to become manifest. The fact is that the 
probability Vs(t) = |a(t)| 2 that an unstable particle <j> 
survives up to time t ~ T is extremely small. Let P^,(t) 
be Pt/>(i)\. T ~ 10~ fe , where k >■ 1, then there is a 
chance to observe some of particles <p survived at t ~ T 
only if there is a source creating these particles in Ms 
number such that Vs(t)\ tr ^ T Ms 3> 1- So if a source 



10' 



exists that creates a flux containing Ms 
stable particles and I ^s> k then the probability the- 
ory states that the number N surv of unstable particles 



-Ps(t)\ t ^ T Ms 
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l-k 



1 , has to survive up 



to time t ~ T. Sources creating such numbers of unstable 
particles are known from cosmology and astrophysics: as 
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example of such a source can be considered processes tak- 
ing place in galactic nuclei (galactic cores), inside stars, 
etc. According to estimations of the luminosity of some 
7-rays sources the energy emitted by these sources can 
even reach a value of order 10 52 [erg/s], 0, Il8l - [20j . and 
it is only a part of the total energy produced there. So, 
if one has a source emitting energy 10 50 [erg/s] then, eg., 
an emission of No — 6.25 x 10 4 ' [1/s] particles of en- 
ergy 10 18 [eV] is energetically allowed. The same source 
can emit J\f ~ 6.25 x 10 56 [1/ s] p articles of energy 10 9 
[eV] and so on. If one follows [13( and assumes that for 
laboratory systems a typical value of the ratio E^/r° is 
E$/r$ > O(10 3 - 10 6 ) and then taking, eg. E%/ T$ = 
10 6 one obtains from Q that A^(T) - 2.5 x 10" 26 A/'o 
and from the estimation of T used in [l3| (see (11), (12) 
in [Hi) th at M^T) ~ Nq. This means that there 

arc AQ,(T) ~ 14 x 10 21 particles per second of energy 
W 4> = 10 i8 [ cV ] or AA (T) - 14 x 10 30 particles of energy 
= 10 9 [eV] in the case of the considered example 
and T calculated using (|7|). On the other hand from 
T obtained for the model considered in [l3| one finds 
W (T) - 6.25 x 10 17 and A^(T) ~ 6.25 x 10 26 respec- 
tively. These estimations show that astrophysical sources 
are able to create such numbers A/"o of unstable particles 
that sufficiently large number N<j,{T) 3> 1 of them has 
to survive up to times T when the effect described above 
should occur. So the numbers of unstable particles pro- 
duced by some astrophysical sources are sufficiently large 
in order that a significant part of them had to survive up 
to the transition times and therefore to emit electromag- 
netic radiation. The expected spectrum of this radiation 
can be very wide: From radio frequencies up to 7-rays 
depending on energy distribution function u(E) of the 
unstable particle emitting this radiation. 

We have shown that charged unstable particles or neu- 
tral unstable particles with non-zero magnetic moment, 
which survived up to transition times or longer, should 
emit electromagnetic radiation. We have also shown 
that only astrophysical processes can generate sufficiently 
huge number of unstable particles in order that this emis- 
sion could occur. From our analysis it seems to be clear 
that the effect described in this letter may have an as- 
trophysical meaning and help explain the controversies, 
which still remain, concerning the mechanisms that gen- 
erates the cosmic microwave, or X-, or 7-rays emission, 
e.g. it could help explain why some space areas (bubbles) 
without visible astronomical objects emit microwave ra- 
diation, X- or 7-rays. Indeed, let us consider active 
galactic nuclei as an example. They emit extremely huge 
numbers of stable and unstable particles including neu- 
trons (see eg. (H) along the axis of rotation of the galaxy. 
The unstable particles, which reached distances d ~ cLt 
from the galactic plane, should emit electromagnetic ra- 
diation. So a distant observer should detect enhanced 
emission of this radiation coming from bubbles with the 



centra located on the axis of the galactic rotation at aver- 
age distances dr from the galactic plane (see Fig. (3)). In 
the case of neutrons dx n can be extremely large. There- 
fore a possible emission of the electromagnetic radiation 
generated by neutrons surviving sufficiently long seems 
to be relatively easy to observe and it should be possible 
to determine dr„ ■ Now having realistic sufficiently accu- 
rate ui(E) for neutrons we are able to calculate T n and 
to find £ n {t) and its local maxima at transition times. 
Thus if the energies W n , (i.e., 72, ), are known then in 
fact we know velocities v n and we can compute dr n and 
distances where £ n (t) has maxima. All these distances 
fix the space areas where the mechanism discussed should 
manifest itself. This suggests how to test this mechanism: 
The computed dx n can be compared with observational 
data and thus one can test if the mechanism described in 
our letter works in astrophysical processes. 

Note that all possible effects discussed in this paper 
are the simple consequence of the fact that the instan- 
taneous energy £${t) of unstable particles becomes large 
for suitably long times compared with £ 9 and for some 
times even extremely large. This property of £ < f>(t) is a 
purely quantum effect resulting from the assumption that 
the energy spectrum is bounded from below and it was 
found by performing an analysis of the properties of the 
quantum mechanical survival probability a(t). 
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